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We investigate the harmonically trapped 2D fermionic systems with a effective spin-orbit coupling 
and intrinsic s-wave superfluidity under the local density approximation, and find that there is a 
critical value for Zeeman field. When the Zeeman field larger than the critical value, the topological 
superfluid phases emerge and coexist with the normal superfluid phase, topological phase separation, 
in the trapped region. Otherwise, the superfluid phase is topologically trivial. 
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1. INTRODUCTION 

Recently, topological supcrconductors/supcrfluids(TSCs/TSs) which cannot be described by the Landau symmetry 
breaking theory, but the topological order [1], have attracted a great attentions in condensed matter physics [2-8]. 
TSCs which resemble topological insulators (TI) [9-12], have gapped bulk and gapless Majorana edge states on the 
boundaries which are protected by the topological properties of the bulk for the bulk-edge correspondence. The distin- 
guishable phases have same symmetry but cannot change into each other without topological phase transition(TPT). 
The realization of TSCs/TSs has been proposed in a large spin-orbit coupling(SOC) semiconductors with a proximity 
induced s-wave pairing potential and ultracold atomic gases with intrinsic superfluidity in the presence of a Zeeman 
field and an effective SOC. 

Ultracold atomic gases with remarkable controllability provide an ideal stage to investigate the phenomena in 
the condensed matter physics [1.3-20]. The effective SOC which is essential for the realization of different TIs and 
TSCs/TSs, can be generated by utilizing the spatial varying laser fields [21-23] and has been realized in experiment 
[24, 25]. The topological properties of the superconductors/supcrfluids were investigated in [5, 26, 27]. The SOC can 
change the ground state of the boson condensate [28, 29], too. 

In this paper, we investigate the topological properties of 2D s-wave supcrfluids trapped in a harmonic potential 
with a SOC in local-density approximation(LDA). The superfluid phase is topologically trivial in the absence of the 
Zeeman field. An effective Zeeman field, breaking the time-reverse symmetry and induced by the spin imbalance of 
the population, is needed to realize the topologically nontrivial phases. The topological nontrivial phases emerge in 
the trapped region when the Zeeman field is larger than the critical value. There is topological phase transition when 
the Zeeman field crosses the critical values. In the presence of the harmonic potential, the chemical potential and the 
pairing gap are functions of the position coordinate. Such that the topological properties of the superfluids are also 
changing with the coordinate. Analogy to the case of without the SOC [15, 20, 30-32], we find there exist topological 
phase scparation(TPS) phenomena which is the coexistence of topological trivial and nontrivial superfluid phases 
in the trapped region. The gapless Majorana edge states are localized at the interfaces between the topologically 
superfluid phases (TS) and the normally superfluid phases (NS) in the trapped region. 

This paper is organized as follows. In Sec. II, first, we introduce the Hamiltonian of the model. Then, we derive the 
energy gap closing conditions and introduce the TKNN numbers. In Sec. Ill, we investigate the topological properties 
of the superfluid phase and find there exist topological phase separation phenomena in the trapped region. We also 
extend the analysis to the continuum model. The conclusions are given in Sec. IV. 
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2. FORMALISM OF THE SYSTEM 

We consider spin-singlet superfluids with the Rashba SOC on a 2D square optical lattice in the harmonic potential, 
which is described by the Hamiltonian(fi, = ks = 1): 

H = H + H S o + H s , (1) 

where Hq is the kinetic term with a Zeeman field, Hso is the spin-orbit interaction term, and H s is the pairing 
potential term. They take the following forms 

H " = - f E E ( c Ui**« + h.c.) - Y^Mr) + h)4 t c, n + ( M (r) - h)^}, (2) 

i,r, a e—x,y i,r 
H SO = -A^[( C Li4 C it-4+ii C it) + i ( c L;)i Ci t" C l+ 2 }| Ci t) + /l - C -]' ( 3 ) 

H s = -j;A(r)(^tk), (4) 

i.r 

where c\ a {ci. a ) denotes the creation (annihilation) operator of the fermion with spin a = at site i = (i x ,i y ), 

A is the strength of SOC and A(r) is the pairing potential of the superfluid. /i(r) = fi — V(r) is the local chemical 
potential, here, fj, is global chemical potential and V(r) = moj 2 r 2 /2 is harmonic potential, h is a Zeeman field. 
Introduce ^ r = {c\. ^, ^, c_k,f, c -k,^) with = e ik l C; , the Hamiltonian, in momentum space, can be rewritten 
as 

H = jE*L^(M*k,r, (5) 
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with 



V -«A(k,r)'cr a -e(k, r) + ft<r» + g k • a* J 

where a = (a x , <r y ) arc Pauli matrices, e(k, r) = — 2t(cos fc x -|-cos k y )—fx(r), and gk = 2A(sin k y , — sin k x ). Diagonalizing 
the Hamiltonian, we find the eigenvalues E(\c, r) of the Hamiltonian satisfy: 



E(k, rf = e(k, rf + h 2 + |g k | 2 + | A(k, r) \ 2 ±2^h 2 (e(k,r) 2 + \ A(k,r) \ 2 ) + e(k, r) 2 |g k | 2 . (7) 

The pairing potential and the global chemical potential can be tuned by the inter-species interaction via the 
Feshbach resonance and the filling number. The pairing potential is a smooth function of the position coordinate. 
We consider the Zeeman field h > 0. In the presence of pairing (A(k, r) ^ 0), the spectrum of our model is gaplcss 
only when |gk| = at four momenta ko = (0, 0), (0,7r), (tt, 0), (w, 7r), and the Zeeman magnetic field must satisfies 
h = v /e(ko,r) 2 + |A.(k ,r)| 2 which 

can be written as: 



hi = A /(4t + M (r))2 + |A(k ,r)|2, k = (0, 0) 

h 2 = y /gg + |A(kp,7)P, k = P[(tt, 0)] , (8) 

h 3 = ^(4i- M (r)) 2 + |A(k ,r)| 2 , k = (tt, tt) 

where P[k] is a set of the vectors obtained from permutations of vector k. 
We know that this BdG Hamiltonian has particle-hole symmetry(PHS) 

U(k,r) - -t x U(-k,r) T t x , (9) 

where t x is Pauli matrix. The classification of time-reversal breaking(TRB) BdG Hamiltonian with PHS is Z class 
[6]. The topological numbers which characterize the topological properties of the superfluid phases are integers. The 
Hamiltonian has four energy bands, only two of them can be gapless and others are always gapped. Only the gapless 
energy bands have contribution to the topological number. We consider the s-wave pairing case with a real pairing 
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gap function A(k, r) — A s (r). then, we can deform the Hamiltonian H into a dual Hamiltonian HP by a unitary 
transformation 

U D (k,r) = DH(k,r)D\ (10) 

where, the dual Hamiltonian and the unitary matrix D is given as: 

H D (\cr\-( A(k,r)-ft<r 2 -ie(k,r)a y -ig k -aa y \ , . 

lK ' rJ {ie(k,r)a v +i S)i aya A(k,r) + ha z )' {1L) 

and 

D= i(} m. ( i2) 

We can adiabatically deform the Hamiltonian by taking e(k, r) = without close the energy gap. Such that the 
dual Hamiltonian is decomposed into the following two 2x2 parts for dual up-spin and down-spin, 



/D rir — ( A s (r) - h 2A(sin k y + i sin k x ) 



11 y 2A(sm k y — i sin fc x ) — A s (r)+ft y 

-uD(V r \-( A s (r) + h 2\(~smk y + ismk x )\ , . 

«uA K > r J ^ 2A(-sinfc y -iainfe) -A 5 (r) - ft ,T 1 j 

When the Zccman field ft > 0, the energy gap of the dual up-spin can vanish while the down-spin is always gapped. 
And when the Zeeman field ft < 0, the reverse case happen. For the simplification, we only consider h > case and 
rewrite the Hamiltonian of the dual up-spin as: 

Wft(k,r) =d(k,r) -a, (15) 

with d(k) = (2Asinfc 2/ , — 2 A sin /e^, A s (r) — h). 

As for the harmonic potential, the topological properties of the supcrfluid phases change with the distance away 
from the trapping center as well as the BdG Hamiltonian. When the harmonic potential is weak the LDA is reliable. 
In order to characterize the topological properties of the phases in the trapped region, the TKNN number Itknn, 
The TKNN number Itknn which is equivalent to the first Chern number, is introduced, 

1 /■„-,„ x ddCk) ddCk) 

7 ™=^r kd «-^ x ^ (i6) 

with d(k) = d(k)/|d(k)|. For the TRB system, the TKNN number is equals to zero for the normal phase and nonzero 
for the topological phase. To simplify the calculation, we take A s (r) = A s . Considering the Hamiltonian with the fixed 
pairing gap A s and a positive tunable Zccman field h, the TKNN number Itknn (ft) can change only if the Hamiltonian 
crosses the gapless point where c? 2 (k, ft) = 0. We find that the phase is topologically trivial(/T.R"Aw(0) = 0) in the 
absence of a Zeeman field. So we only need to calculate the change in Itknn (ft) at critical points which arc ft = A s at 
the momenta k c = (0, 0), P\(if, 0)], (it, it). We must keep in mind that the different momenta k c indicate the different 
adiabatic deformation e(k c , r) = 0. Such that the three critical momenta correspond to three gapless conditions. 

Around the critical momentum k c = (0, 0), the vector d(k) has an approximate form d(k) = (2\k y , —2Xk x , A s (r) — 
ft). Taking a cutoff A < tt and in the limit ft — A s = Sh -C A s , the TKNN number expression can be divided into low 
and high energy parts, 

Itknn * Afl^ • ^ x ^ 



4tt J d 3 (k) dk x dk y 
1 , f „ f „ . A 2 cos k x cos k v 5h 

-(/ dk+ I dk) — ; \ - * 



K J\k\<A J\k\>\ {4\ 2 {sin 2 k x + sin 2 k y ) + Sh 2 ) 2 
= 4™(^A) + 4™(^A). (17) 

Since there is no gapless point in the high momentum region |k| > A, the change in Itknn only come from the low 
energy part Itknn- The change in Itknn can be calculated at critical momentum k c = (0, 0), 

AA1 TKNN\Sh=Q- — x > \ la ) 
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FIG. 1: We show the phase diagrams on the h — r plane for the pairing gap A s = t, and the global chemical potential (a) 
/i = — 2t , (b) /i = 0, and (c) fj,(r) = 2t. The phases of the shadow regions are TS and others are NS. The TKNN number 
Itknn is ±1 for the slight shadow regions and 2 for the dark shadow regions. The entire trapped region is NS phase for 
h/t < 1. As the Zeeman field increasing, there is topological phase transition at h/t — 1. The TS phases emerge and coexist 
with the NS for the Zeeman field h/t > 1. 

Similarly, the change in TKNN number at other critical points can be obtained and given as: 

f -1, k c = (0,0) 
AI TKNN \ s s 'Ct- = \ 2, k c = P[(0 j7 r)] , (19) 
[ -1, k c = (tTjTt) 

3. TOPOLOGICAL PHASE SEPARATION 

In the presence of the harmonically trapping potential, the pairing gap is a function of the distance away from the 
trapped center which can be self-consistently determined with the global chemical potential via the minimization of 
the free energy. By using the change in Itknn at the critical points, we can study the topological properties of the 
supcrfiuid phase in the trapped region. The change in Itknn at the boundaries can be listed as below: 

f -1, h = hi 

AlTKNN\ S s h h Zo- = \ 2 > h = h 2 . (20) 
[-1, h = h 3 

In the absence of the Zeeman field, the TKNN number vanish such that the superfluid phases is topological trivial 
in the entire trapped region. As the Zeeman field increasing, the TKNN number change to nonzero at the critical 
Zeeman field h c = min(/ii, h%, /13). Such that the topologically nontrivial phase emerge in the trapped region. There 
is topological phase transition when the Zeeman field crosses the critical value. When the Zeeman field is larger than 
the critical value, the topologically trivial and nontrivial phases coexist in the trapped region. This topological phase 
separation is analogy to the phase separation phenomena [32] in the absence of SOC case. 

The interfaces between the superfluid phases of different topological properties are localized in the trapped region. 
From the bulk-edge correspondence, a topologically nontrivial bulk guarantees the existence of topologically stable 
gapless edge states on the interfaces. For this 2D superfluids breaking the time-reversal symmetry but having the 
PHS, the gapless edge states are chiral Majorana fermion modes. The number of gapless edge states is defined by 
the TKNN number Itknn- For the odd Itknn phase, there are odd Majorana fermions in each vortices which are 
neither bosons nor fermions but non-Abelian anyons and obey the non-Abelian statistics. For the even Itknn phase, 
as for the vortices having even Majorana fermions, there are no non-Abelian anyons. 

Now, we assume the pairing gap is independent of the distance away from the trapped center (A s (r) = A s ) and the 
harmonic potential is V(r) = mw 2 r 2 /2 = O.ltr 2 where r is in the unit of the lattice space. We can obtain the phase 
diagrams on the h — r plane with the s-wave pairing gap A s = t for the different global chemical potential, shown in 
Fig.[l]. 

From Fig.[l], we find that the Zeeman field has a critical value at h = t. When the Zeeman field h < t, the 
superfluid phase is topological trivial in the entire trapped region. As the Zeeman field increasing, there is TPT at 
h/t = 1. When the Zeeman field h > t, the topologically superfluid phases emerge and coexist with the normally 
superfluid phases in the trapped region. The superfluid phases are topologically nontrivial for the shadow regions 



(a) (b) 

FIG. 2: We show the diagrams of different types TPS phenomena in real space. The TPS of (a) is N-T-N and (b) is T-N case 
in the trapped region. The phases in the shadow regions are the TS, and the other regions are the TS. 
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FIG. 3: We show the phase diagrams on the h — r plane for the continuum model with the pairing gap A s = O.IEf, and the 
global chemical potential (a) /i = — 0.3-Ef and (b) /i = O.SEf- The phases of the shadow regions are TS and others are NS. 
The entire trapped region is NS phase for h/t < 1. The TS phase emerges and coexists with the NS for the Zeeman field 
h > minrWfi(r) 2 + \A s (r)\ 2 ). 

with the others are topologically trivial. The TKNN number Itknn is ±1 for the slight shadow regions and 2 for 
the dark shadow regions. When the global chemical potential /i = — 2t, shown in Fig. [la], there is TPS phenomena in 
the trapped region for h > t. The TS circle is immersed in the NS disk for the Zeeman field 1 < h/t < 2 and h/t > 6. 
The phase diagram in the real space is N-T-N as shown in Fig. [2a]. When the Zeeman field 2 < h/t < 6, the TS disk 
is surrounded by NS circle shown in Fig. [2b]. For /j, = case, there is TS disk with Itknn = 2 when the Zeeman 
field 1 < h/t < 4, shown in Fig. [lb], and the TPS phenomena occur only when the Zeeman field h > t. Fig.[lc] gives 
the phases diagram on the h — r plane for the global chemical potential = 2t. The TS phases emerge and coexist 
with the NS phase only when h > t, too. 

In Fig. [2], we show two types TPS phenomena in the real space. There are two interfaces between the NS and TS 
in Fig. [2a] and one in Fig. [2b]. The interfaces are Majorana gaplcss edge states. 

The TPS is a ubiquitous phenomena in the weak trapped system with SOC. We extend the above discussion 
to the continuous model. For the continuum s-wave case, the Hamiltonian can be written in the form of Eq.(5) 
with e(k, r) = ti 2 k 2 /2m — n(r), gk = a(k y ,—k x ) and A(k, r) = A s (r). The gapless conditions are k = and 
h = y/v(r) 2 + |A s (r)| 2 . There is TPT at the critical Zeeman field (h c = min r (^/ ^(r) 2 + |A s (r)| 2 )). The TS emerge 
and coexist with the NS when h > h c . The topological phase diagrams on h — r plane are shown in Fig. [3a] for 
fi = -0.3E F and Fig. [3b] for fi = 0.3E F . 

4. CONCLUSIONS 

In this paper, we investigate the trapped 2D fcrmionic systems with a effective SOC and intrinsic s-wave superfluid- 
ity, and find that the Zeeman field exists a critical value. There is topological phase transition at the critical Zeeman 
field. When the Zeeman field larger than the critical value, the topological superfluid phase emerge and coexist 
with the normal superfluid phase in the trapped region which we called topological phase separation. Otherwise, the 
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superfluid phase is topologically trivial. We extend the analysis to the continuum case and find that the Zeeman field 
also exists a critical value. For the continuum case, the topologically superfluid phases emerge and coexist with the 
normally superfluid phase When the Zeeman field larger than the critical value. Otherwise, the superfluid phase is 
topologically trivial, too. 
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